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TORSIONAL  SHEAR  STRENGTH  OF  FILAMENT- WOUND 
GLASS-EPOXY  TUBES 

By  L.  David  Wall,  Jr.,  and  Michael  F.  Card 
Langley  Research  Center 

SUMMARY 

Results  are  presented  from  torsion  tests  conducted  on  36  multilayered,  filament- 
wound,  glass-epoxy  tubes.  Configurations  with  helical  windings  and  with  alternating  heli¬ 
cal  and  circumferential  windings  were  investigated  for  various  winding  angles.  Under 
small  loadings,  shear  moduli  deduced  from  linear  shear  stress-strain  curves  were  found 
to  be  in  reasonable  agreement  with  analytical  predictions.  Under  larger  loadings,  various 
degrees  of  nonlinearity  in  shear  stress-strain  curves  were  encountered,  depending  on  the 
helical  winding  angle.  Experimental  torsional  strengths  were  defined  by  a  0.2-percent 
offset  yield  stress  or  by  maximum  stress  when  large  nonlinearities  did  not  exist.  These 
strengths  were  compared  with  torsional  buckling  predictions  for  orthotropic  cylinders  and 
with  material  strength  predictions  based  on  orthotropic  yield  criteria  and  elastic  stress 
analysis.  Computed  elastic  buckling  stresses  were  considerably  higher  than  the  experi¬ 
mental  strengths  for  most  of  the  test  specimens  except  for  those  with  only  30°  and  45° 
windings.  Experimental  torsional  strengths  were  found  to  correlate  with  conventional 
yield  predictions  if  predicted  yielding  in  certain  layers  were  ignored,  or  if  unrealistically 
large  transverse  tensile  and  shear  strengths  of  unidirectional  laminae  were  employed  in 
the  analysis. 


INTRODUCTION 

The  development  of  high-performance  glass,  boron,  and  carbon  filaments  has 
resulted  in  many  studies  of  the  application  of  fiber-reinforced  structures  to  aerospace 
structures.  Successful  design  with  such  structures  requires  a  thorough  knowledge  of  the 
strength  of  fibrous  composites  for  a  wide  variety  of  configurations  under  tensile,  com¬ 
pressive,  and  shear  loads.  Although  some  experimental  information  is  available  on  the 
strength  of  fibrous  composites,  there  is  no  substantial  body  of  reliable  data  available  for 
shear  loadings. 

In  the  experimental  determination  of  the  shear  strength  of  a  multilayered  composite, 
a  variety  of  test  methods  including  beam,  plate,  and  tube  tests  can  be  employed.  The 
circular  tube  or  cylindrical  shell  is  considered  to  be  the  most  reliable  specimen  for 


measurement  of  shear  strength  since  torsional  loadings  produce  the  most  easily  under¬ 
stood  state  of  stress  in  the  test  section.  Unfortunately,  test  programs  on  the  torsional 
shear  strength  of  tubular  specimens  have  been  generally  limited  to  a  few  isolated  tests. 
Available  test  data  for  glass-  and  boron- reinforced  tubes  or  cylindrical  shells  are  given 
in  references  1  to  7. 

The  purpose  of  the  present  paper  is  to  report  the  results  of  an  experimental  investi¬ 
gation  of  the  torsional  strength  and  stiffness  of  36  glass-epoxy  multilayered  filament- 
wound  tubes.  Two  basic  winding  configurations  were  studied:  one  with  helical  windings, 
the  other  with  alternating  helical  and  circumferential  windings.  Torsional  stiffnesses  and 
strengths  of  tubes  with  various  helical  winding  angles  were  determined  experimentally. 
Shear  moduli  are  compared  with  analytical  estimates.  Shear  strengths  based  on  ultimate 
or  yield  stress  are  compared  with  predictions  from  orthotropic  buckling  theory  and  with 
orthotropic  yield  criteria  using  the  methods  summarized  in  reference  8.  The  experimen¬ 
tal  program  described  herein  has  been  complemented  by  a  similar  program  (unpublished) 
in  which  the  axial  compressive  strength  of  glass-epoxy  tubes  was  investigated. 

SYMBOLS 

Measurements  and  calculations  were  made  in  U.S.  Customary  Units.  The  values  for 
the  physical  quantities  are  given  both  in  the  International  System  of  Units  (SI)  and  in  U.S. 
Customary  Units.  A  table  of  conversion  factors  is  given  in  appendix  A.  The  relationship 
between  these  two  systems  of  units  can  be  found  in  reference  9. 

a,b  constants  associated  with  orthotropic  yield  criteria  (eq.  (5)) 

Bij  elastic  constants  in  inverted  form  of  strain-stress  relations  (B3) 

Cij  extensional  stiffness  of  cylinder  wall 

Dy  bending  stiffness  of  cylinder  wall 

E  Young’s  modulus 

Young’s  modulus  of  unidirectional  layer  in  direction  parallel  and 
perpendicular  to  fibers,  respectively 

Ex,Ey  Young's  moduli  for  layer  associated  with  cylinder  surface  coordinates 
Ex  effective  axial  modulus  of  cylinder  wall 
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G 

gLT 

Gxy 

Gxy 

Kij 

l 

M  x 

mi,m2 

Nx,Ny 

Nxy 

Nxy 

P 

R 

S 

T 

t 

u,v,w 

V 


shear  modulus 

shear  modulus  of  unidirectional  layer 

shear  modulus  for  layer  associated  with  cylinder  surface  coordinates 
effective  shear  modulus  of  cylinder  wall 

stiffness  of  cylinder  wall  associated  with  bending- extension  coupling 
length  of  cylinder 
bending- moment  resultant 

auxiliary  constants  associated  with  shear  stress  in  a  generally  orthotropic 
layer  (see  eq.  (Bl)) 

axial  and  circumferential  load  resultants 
shear  stress  resultant 

applied  shear  stress  resultant  at  boundary  of  shell 
internal  pressure 
radius  of  cylinder 

allowable  yield  strength  in  shear  for  unidirectional  layer 
applied  torque 
thickness  of  tube  wall 

displacements  of  cylinder  wall  associated  with  x,  y,  and  z  coordinates, 
respectively 

volume  fraction  (always  subscripted),  ratio  of  constituent  volume  to 
total  volume 
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allowable  yield  strengths  in  layer  parallel  and  transverse  to  fibers, 
respectively 

surface  coordinates,  measured  parallel  to  cylinder  axis,  perpendicular 
to  cylinder  axis  and  normal  to  cylinder  surface,  respectively  (see  fig.  1) 

orthotropic  curvature  parameter 

helical  winding  angle  measured  from  axis  of  cylinder  (see  fig.  1) 
shearing  strain 
normal  strain 

correlation  factors  associated  with  shear  and  extension,  respectively 
(see  eqs.  (1)) 

change  in  curvature 

Poisson's  ratio 

Poisson's  ratios  of  a  layer  associated  with  extension  induced  by  inplane 
loading  parallel  and  perpendicular  to  fiber  direction,  respectively 

Poisson's  ratios  of  a  layer  associated  with  extension  induced  by  loading 
parallel  and  perpendicular  to  axis  of  cylinder,  respectively 

effective  Poisson's  ratio  of  cylinder  wall  associated  with  axial  loading 

density  of  material 

normal  stress 

shear  stress 


applied  average  shear  stress 


Subscripts: 


B  quantities  associated  with  buckling 

f  fiber  property 

i  ith  layer 

m  matrix  property 

L  fiber  direction 

T  direction  transverse  to  fiber 

v  void 

max  maximum 

x  x-direction 

y  y-direction,  yield 

Superscript: 

T  total 

A  subscript  preceded  by  a  comma  denotes  differentiation  with  respect  to  the  sub¬ 
script.  Bars  over  symbols  denote  applied  loads  at  the  boundary. 

EXPERIMENTAL  INVESTIGATION 

Test  Specimens 

The  winding  geometry  of  the  test  specimens  of  the  present  investigation  is  indi¬ 
cated  in  figures  1  and  2.  Figure  1  shows  the  convention  used  in  describing  the  winding 
angles  (±or)  as  well  as  the  basic  coordinate  system  (x  and  y)  and  the  sign  convention 
for  a  positive  torque  T.  The  two  types  of  winding  configurations  investigated  are  shown 
in  the  photomicrographs  presented  in  figure  2.  The  line  of  sight  for  the  figure  was  the 
y-axis  so  that  cross  sections  of  circumferentially  wound  filaments  project  as  circles. 

The  configuration  with  all  helical  wraps  (fig.  2(a))  is  composed  of  10  layers  wrapped  in  the 
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sequence  +a,  -a,  +a,  -a,  +a,  -a,  +a,  -a,  +«,  -a.  The  configuration  with  alternating 
helical  and  circumferential  wraps  (fig.  2(b))  is  composed  of  12  layers  wound  in  the 
sequence  +a,  -a,  90,  +<*,  -a,  90,  +a,  -a,  90,  +a,  -a,  90.  In  both  configurations,  all 
layers  were  approximately  equal  in  thickness.  An  unusual  feature  of  the  test  specimens 
was  that  they  did  not  contain  filament  crossovers  or  weaving  which  customarily  occur  in 
the  filament  winding  of  cylinders  (see,  for  example,  ref.  10,  p.  50);  thus,  each  layer  had 
only  unidirectional  filaments.  The  two  photomicrographs  shown  in  figure  2  also  show  the 
irregular  fiber  spacing  and  the  white,  resin- rich  areas  which  are  typical  of  glass  filament- 
wound  structures.  The  dark  areas  in  the  figure  are  believed  to  have  resulted  from  nonuni¬ 
form  polishing  of  the  specimen  and  may  be  an  indication  of  void  content. 

The  test  specimens  were  wet  wound  from  E  glass  and  ERL  2256  epoxy  as  long  tubes 
which  were  cured  at  elevated  temperature.  The  properties  of  the  constituent  materials 
and  the  cure  cycle  are  given  in  table  I.  The  tabulated  values  of  matrix  modulus  and 
matrix  density  were  obtained  from  7.6-cm  (3-inch)  diameter  by  15.2-cm  (6-inch)  long 
blocks  of  resin  which  were  tested  in  compression.  The  density  and  modulus  of  glass  and 
Poisson's  ratios  shown  are  representative  values  for  these  materials.  Four  test  speci¬ 
mens  were  cut  from  each  long  tube.  Two  of  the  specimens  were  employed  in  the  investi¬ 
gation  described  herein.  The  other  two  specimens  were  used  in  a  compressive  strength 
investigation.  Each  test  specimen  had  the  following  nominal  dimensions:  7.6-cm  (3-inch) 
inside  diameter,  18- cm  (7-inch)  total  length,  and  a  total  wall  thickness  of  0.152  cm 
(0.06  inch).  There  were  24  helically  wrapped  specimens  with  wrap  angles  varying  from 
15°  to  90°  in  15°  increments  and  12  specimens  with  alternating  helical  and  circumferen¬ 
tial  wraps  with  wrap  angles  of  15°,  30°,  and  45°. 

Measured  geometry  and  constituent  volume  fractions  for  each  test  specimen  are  pre¬ 
sented  in  tables  II  and  III.  The  two  torsion  specimens  cut  from  each  long  tube  are  num¬ 
bered  successively  as  1  and  2,  3  and  4,  and  so  forth  in  the  tables.  The  helical  angles  ±a 
and  the  circumferential  angles  on  each  specimen  were  measured  from  photographs 
obtained  by  wrapping  photographic  film  around  the  circumference  of  the  long  tube  and  illu¬ 
minating  the  inside  of  the  tube.  The  values  of  wall  thickness  t  shown  are  the  average  of 
several  random  measurements  on  each  specimen.  Scatter  in  thickness  measurements  was 
generally  about  ±5  percent  of  values  listed  in  tables  II  or  III  although,  in  a  few  specimens, 
deviation  of  ±6  to  ±12  percent  existed.  (See  values  of  t  marked  with  an  asterisk  in 
tables.) 

The  tabulated  values  of  fiber  volume  fraction  Vf  and  matrix  volume  fraction  vm 
are  averages  obtained  for  four  coupons  cut  from  the  walls  of  each  test  specimen.  In  order 
to  determine  Vf  and  vm  in  a  manner  such  that  void  volume  fraction  could  be  estimated, 
the  total  volume  of  each  coupon  was  obtained  by  comparing  the  weight  of  the  coupon  in  air 
with  the  weight  of  the  coupon  submerged  in  water.  Then  burnout  tests  of  3-hour  duration 
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at  866  K  (1100°  F)  were  conducted  on  each  coupon  to  obtain  its  glass  and  matrix  material 
content.  The  difference  between  unity  and  the  sum  of  average  glass  and  average  matrix 
volume  fractions  was  used  as  an  estimate  of  void  content.  Scatter  in  Vf  and  vm 
between  coupons  of  the  same  specimen  was  about  ±3  percent  of  the  values  in  tables  II  and 
HI.  Since  the  void  content  vv  was  obtained  as  the  difference  of  two  large  numbers,  the 
tabulated  values  are  probably  accurate  only  to  the  nearest  hundredth. 

Test  Procedures 

To  prepare  the  specimens  for  torsion  testing,  the  ends  of  each  specimen  were  rein¬ 
forced  with  a  0.17-  by  5-cm  (1/16-  by  2-inch)  aluminum -alloy  split  ring  bonded  to  the 
external  surface  of  the  tube  with  a  room-temperature  curing  epoxy.  The  addition  of  end 
reinforcement  resulted  in  a  test-section  length  of  about  7.6  cm  (3  inches).  The  specimens 
were  tested  in  a  torsion  testing  machine  of  6.8  kN-m  (60  000  in- lb)  capacity  at  the  Langley 
Research  Center.  One  end  of  each  specimen  was  mounted  to  the  stationary  head  of  the 
machine.  The  other  end  was  attached  to  the  rotating  head  which  induced  a  torque  deter¬ 
mined  by  a  sensitive  weighing  system.  Before  each  test,  the  gripping  chucks  of  the  tor¬ 
sion  machine  were  carefully  alined  to  ensure  that  the  loading  transferred  to  the  test  speci¬ 
men  would  be  pure  torsion.  The  specimen  was  adapted  to  the  chucks  by  inserting  an 
accurately  machined  plug  in  each  end  to  a  depth  of  about  4.8  cm  ^1^-  inches^.  The  plugs 
were  carefully  shimmed  for  axial  alinement  and  were  attached  to  the  specimen  by  screws. 
The  aluminum  inserts  and  split  rings  can  be  seen  in  the  test  setup  shown  in  figure  3. 

To  measure  shearing  strains,  four  wire-resistance-type  strain  rosettes  were  bonded 
to  the  external  surface  of  the  tube.  Three-gage  45°  rosettes  were  employed  and  the  cen¬ 
tral  gage  was  alined  parallel  to  the  axis  of  the  tube.  (See  fig.  3(b).)  The  two  orthogonal 
gages  were  used  to  measure  the  principal  normal  strains  whereas  the  central  gage  was 
used  to  detect  any  axial  strains.  In  order  to  determine  the  relative  angle  of  twist  of  two 
sections  normal  to  the  tube  axis,  most  of  the  specimens  were  equipped  with  two  parallel 
steel  rods.  Each  rod  (0.48-cm  (3/16-inch)  diameter)  was  mounted  in  two  tightly  fitting 
holes  drilled  through  the  walls  of  the  specimen  in  such  a  manner  that  the  mounted  rods 
passed  diametrically  through  the  specimen.  The  two  rods  were  spaced  about  6.4  cm 
| '2 i  inches)  apart.  The  relative  angle  of  twist  of  one  rod  with  respect  to  the  other  was 
determined  by  measuring  deflections  of  the  ends  of  the  rods  with  four  deflection  trans¬ 
ducers  (linear  variable  differential  transformers).  In  figure  3(a)  the  two  transducers  on 
the  near  side  of  the  specimen  are  shown. 

After  completion  of  tests  of  the  45°  tubes  with  circumferential  wraps,  a  review  of 
the  data  suggested  that  some  of  the  specimens  might  be  developing  buckles  before  failure. 
Therefore,  eight  of  the  tubes  (one  from  each  of  the  remaining  configurations)  were  instru¬ 
mented  with  20  additional  single-element  strain  gages.  These  gages  were  located  in 
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back-to-back  pairs  on  the  inner  and  outer  surfaces  of  the  tube  and  were  alined  with  the 
direction  of  principal  compressive  stress  so  that  they  would  act  as  buckle  detectors.  In 
order  to  obtain  strengths  of  the  tubes  in  the  absence  of  stress  concentration  effects,  these 
eight  tubes  were  not  drilled  to  accommodate  the  rods  for  measurements  of  the  relative 
angle  of  twist. 

Two  types  of  loading  tests  were  conducted.  Because  some  nonlinearity  in  shear 
stress-strain  behavior  was  anticipated,  a  preliminary  test  up  to  a  maximum  torque  of 
68  N-m  (600  in-lb)  was  conducted  on  each  specimen  to  determine  the  initial  shearing 
modulus.  Then,  a  test  to  failure  was  conducted  with  a  single  cycle  of  load.  The  loading 
rate  for  each  test  was  essentially  constant  except  at  the  onset  of  initial  loading  and  during 
failure.  These  rates  were  approximately  0.11  N-m/s  (60  in-lb/min)  for  the  preliminary 
test  and  1.1  N-m/s  (600  in-lb/min)  for  the  test  to  failure.  During  the  tests,  data  were 
recorded  at  a  virtually  continuous  rate  on  the  Langley  central  digital  data  recording 
system. 


Test  Results 

Stress-strain  behavior.-  The  data  obtained  during  the  tests  were  reduced  and  plotted 
by  computer  as  shear  stress-strain  curves  for  each  specimen.  A  sample  plot  of  strain 
data  taken  from  the  orthogonal  gages  of  the  rosettes  during  the  failure  tests  of  specimen 
number  2  is  presented  in  figure  4.  The  magnitude  of  the  applied  shear  stress  Txy  is 
plotted  against  the  shear  strain  yxy*  The  applied  shear  stress  Txy  was  obtained  from 
the  applied  torque  r,  by  using  Bredt’s  formula  with  the  tube  mean  radius,  and  the  shear 
strain  yx y  was  obtained  by  doubling  the  magnitude  of  the  orthogonal  strain- gage  reading. 
Data  from  the  central  axial  gage  of  the  rosettes  indicated  that  for  all  practical  purposes, 
a  state  of  pure  shear  was  induced  in  the  test  section  of  each  specimen.  In  figure  4,  the 
curves  depict  data  from  the  orthogonal  gages  plotted  so  that  tensile  strain  gages  alternate 
with  the  compression  gages;  each  pair  of  gages  is  from  the  same  rosette.  Note  that  the 
gages  have  a  nonlinear  response  which  suggests  that  some  yielding  (or  degradation  of  the 
composite  material)  has  taken  place.  Various  degrees  of  nonlinearity  in  such  stress- 
strain  curves  were  observed,  depending  on  the  helical  wrap  angle  of  the  specimen. 

A  further  illustration  of  the  relative  nonlinear  behavior  of  the  test  specimens  is 
given  by  the  shear  stress-strain  curves  presented  in  figure  5.  Those  curves  are  plots 
showing  shear  strain  obtained  from  the  data  of  a  typical  single  orthogonal  gage  of  a  rosette 
for  one  specimen  in  each  a-group.  For  the  helically  wound  specimens,  considerable  non¬ 
linearity  occurred  in  the  15°,  75°,  and  90°  test  specimens  as  illustrated  in  figure  5(a). 

For  the  helically  and  circumferentially  wound  specimens,  the  15°  tubes  showed  most  of 
the  nonlinearity.  (See  fig.  5(b).)  The  major  difference  between  the  two  winding  configura¬ 
tions,  however,  is  the  maximum  shearing  strain  measured  in  configurations  with  the  same 
helical  winding  angle.  Note  from  figures  5(a)  and  5(b)  that  for  specimens  with  15°  and  30° 
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helical  windings,  those  with  circumferential  windings  were  able  to  sustain  about  twice  the 
maximum  shearing  strain  experienced  by  the  all  helically  wound  specimens.  Thus,  the 
ability  of  helically  wrapped  tubes  to  withstand  large  nonlinear  strains  under  torsional 
loading  appears  to  be  considerably  enhanced  by  the  addition  of  interspersed  circumferen¬ 
tial  wraps. 

Similar  trends  in  stress-strain  behavior  were  observed  by  studying  the  relative  angle 
of  twist  as  determined  from  the  deflection  transducer  data.  The  correlation  of  strain- gage 
data  with  reduced  deflection  transducer  data  indicated  that  the  strain- gage  data  were  a 
good  overall  measure  of  specimen  behavior.  Because  of  rigid-body  translations  and  rota¬ 
tions  resulting  from  apparent  slippage  and  testing  machine  motor  vibrations,  the  response 
of  individual  deflection  transducers  was  not  well  behaved.  The  rigid-body  deflections  of 
the  steel  rods  were  on  the  order  of  4  to  5  times  the  relative  deflection  and  thus  the  rela¬ 
tive  rotation  measurements  represent  small  differences  of  large  numbers.  It  is  felt 
therefore  that  the  relative  rotation  data  are  less  accurate  than  the  strain-gage  data,  and 
hence  are  not  presented. 

The  average  initial  slope  of  stress- strain  curves  of  the  type  shown  in  figure  4  for 
failure  tests  were  compared  with  similar  plots  from  the  preliminary  tests  at  low  loadings. 
The  preliminary  test  data  were  linear  and  the  stress-strain  curve  slopes  were  in  good 
agreement  with  the  failure  test  data  as  well  as  with  the  relative  rotation  data.  Average 
values  of  the  initial  slope  of  shear  modulus  Gxy  are  given  in  tables  IV  and  V  for  the 
two  winding  configurations.  Strain  values  employed  in  computing  modulus  were  corrected 
for  the  transverse  gage  sensitivity  of  the  rosette  by  using  recommended  factors  supplied 
with  the  gages.  Scatter  in  modulus  indicated  by  individual  gages  was  about  ±8  percent  of 
the  values  shown  in  the  tables.  Also  shown  in  tables  IV  and  V  are  the  average  applied 
shear  stress  Ty  corresponding  to  yielding  as  conventionally  defined  in  uniaxial  tension 
tests  of  metals.  The  values  were  obtained  by  constructing  the  linear  stress-strain  curve 
offset  at  0.2-percent  shear  strain  on  each  of  the  stress-strain  curves  (see  fig.  4)  and 
averaging  the  values  of  stress  at  the  intersections  with  the  strain-gage  data  curves.  The 
absence  of  values  for  Ty  in  tables  IV  and  V  indicates  that  strain  values  for  the  specimen 
at  failure  were  not  large  enough  to  define  the  offset  yield. 

Failure.-  The  test  specimens  failed  catastrophically  at  the  maximum  torques  Tmax 
and  shear  stresses  Tmax  presented  in  tables  IV  and  V.  Two  types  of  fracture  patterns 
were  observed  on  the  test  specimens  as  shown  in  figure  6.  In  the  helically  wound  speci¬ 
mens,  a  line  of  separation  occurred  that  ran  parallel  to  the  filaments  loaded  in  tension 
(the  +a  line  of  fig.  1).  The  fracture  line  can  be  seen  on  the  tubes  with  60°  and  75° 
windings  shown  in  figure  6.  For  the  tubes  with  circumferential  windings,  the  fracture 
pattern  was  a  band  rather  than  a  line  and  complete  separation  was  not  evident.  The  direc¬ 
tion  of  the  band  was  only  approximately  parallel  to  the  helical  filaments  stressed  in 
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tension.  The  fracture  band  can  be  seen  in  the  tubes  with  15°  windings  shown  in  figure  6  by 
observing  the  disrupted  lines  of  reflected  light  and  the  areas  of  lighter  color.  The  lighter 
color  gives  evidence  of  delamination  occurring  under  the  external  surface  of  the  tube.  A 
close-up  view  of  the  appearance  of  the  failure  band  under  loading  is  given  in  figure  3(b). 
The  band  appears  as  a  ripple  line  so  that  filaments  with  negative  winding  angles  {-a  in 
fig.  1)  are  wrinkled.  In  figure  3(b)  failure  appears  to  emanate  from  one  of  the  holes 
accommodating  a  steel  rod.  However,  in  many  other  instances,  the  failure  line  or  band  did 
not  include  the  rod  holes.  Furthermore,  the  specimen  without  holes  in  the  same  a-group 
had  a  shear  strength  similar  to  specimens  with  holes.  Thus,  there  appears  to  be  no  exper¬ 
imental  evidence  that  stress  concentration  effects  introduced  by  the  rod  holes  had  any  del¬ 
eterious  effect  on  the  torsional  strength  of  the  tubes. 

In  some  of  the  specimens,  a  buckle  pattern  was  observed  before  the  final  fracture 
pattern  developed.  Additional  evidence  of  buckling  was  obtained  by  observing  strain 
reversal  in  the  strain  gage  data  from  the  specimens  which  were  instrumented  to  detect 
buckling  deformations.  An  example  of  stress-strain  data  with  strain-reversal  trends  is 
presented  in  figure  7.  The  data  shown  are  for  specimen  number  10  (a  =  45°)  and  are 
taken  from  strain  rosettes  plotted  in  the  same  manner  as  in  figure  4.  From  a  comparison 
of  the  data  of  figures  4  and  7,  it  appears  that  buckling  occurred  in  specimen  10  whereas 
specimen  2  experienced  a  material  failure.  On  the  basis  of  studies  of  stress-strain 
behavior  near  failure,  specimens  in  which  significant  strain- reversal  behavior  occurred 
are  identified  in  tables  IV  and  V. 

An  additional  observation  during  the  loading  history  of  the  test  specimens  should  be 
noted.  Each  specimen  emitted  several  cracking  sounds  during  the  strength  test.  In  some 
cases  these  noises  were  accompanied  by  a  drop  in  applied  load  which  upon  further  twisting 
was  followed  by  total  recovery  of  load  along  the  same  stress- strain  curve.  For  a  few 
tests  the  noise  level  and  specimen  acceleration  were  recorded;  however,  there  was  no  cor¬ 
relation  between  these  measurements  and  the  stress-strain  characteristics  already  dis¬ 
cussed.  Other  investigators  have  reported  similar  unsuccessful  correlation  attempts.  It 
seems  that,  at  present,  there  is  no  satisfactory  method  of  assessing  the  significance  of 
this  localized  cracking  on  the  structural  integrity  of  composite  materials. 

COMPARISON  OF  TEST  DATA  WITH  ANALYSIS 
Structural  Stiffness 

The  structural  stiffnesses  of  the  tubes  were  obtained  analytically  from  an  analysis  of 
the  stresses  and  deformations  in  a  finite- length  generally  orthotropic  cylindrical  shell 
loaded  in  torsion.  The  term  "generally  orthotropic"  is  used  to  describe  the  character  of 
the  helical  layers  at  +a  and  -a  which  are  such  that  their  principal  axes  of  orthotropy 
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(the  directions  parallel  and  transverse  to  the  fibers)  are  not  parallel  or  perpendicular  to 
the  cylinder  axis.  (See,  for  example,  ref.  1.)  The  analysis  revealed  that  the  effects  of 
restraints  at  the  boundary  of  the  shell  should  be  very  localized;  therefore,  a  membrane 
state  of  stress  should  yield  accurate  stresses  and  deformations  over  most  of  the  shell. 

The  details  of  the  analysis  are  outlined  in  appendix  B. 

In  order  to  estimate  the  structural  stiffnesses  of  the  tubes  from  the  properties  of  the 
glass  and  epoxy  constituents  (see  table  I),  the  elastic  constants  along  coordinates  associ¬ 
ated  with  the  filaments  were  computed  from  the  Halpin-Tsai  equations  discussed  in  refer¬ 
ence  8.  These  equations  are  relatively  simple  to  employ  and  have  been  empirically  cor¬ 
related  with  trends  from  elasticity  solutions  treating  idealized  arrays  of  inclusions  in  a 
matrix  (for  example,  ref.  1).  If  the  subscript  L  denotes  a  direction  parallel  to  the  fiber 
axes  and  the  subscript  T,  a  direction  perpendicular  to  the  fibers,  the  equations  employed 
for  elastic  constants  for  a  unidirectional  layer  can  be  written  as 


with 
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In  equations  (1),  tq  and  ?s  are  empirical  correlation  factors  associated  with 
transverse  normal  load  and  shear  load  respectively;  the  values  assigned  to  these  param¬ 
eters  were  selected  to  agree  with  experimental  results  of  the  present  investigation. 

As  outlined  in  appendix  B,  the  elastic  constants  for  the  unidirectional  layer  are 
transformed  into  elastic  constants  associated  with  the  x,y  coordinates.  (See  fig.  1.) 

The  structural  stiffnesses  of  the  shell  are  then  obtained  by  inversion  of  the  strain-stress 
law  and  integration  over  the  total  tube  wall  thickness.  The  effective  shear  moduli  of  the 
wall  Gxy  of  the  tube  can  be  computed  from  the  integrated  shear  stiffness.  (See 
eq.  (B19).)  A  comparison  of  the  analytically  and  experimentally  determined  values  of  the 
shear  modulus  is  presented  in  figure  8.  The  calculated  curves  were  computed  for  a  nomi¬ 
nal  glass  volume  fraction  of  0.63  and  are  compared  with  data  with  similar  volume  frac¬ 
tions.  Calculations  for  specific  specimens  with  actual  volume  fractions  are  presented  in 
tables  IV  and  V.  Based  on  the  results  shown,  it  can  be  seen  that  the  computed  results 
are  in  good  agreement  with  experimental  values,  the  maximum  deviation  being  about 
10  percent. 


Torsional  Buckling 

As  previously  mentioned,  certain  experimental  observations  suggested  that  some  of 
the  test  specimens  may  have  buckled  prior  to  final  fracture.  The  torsional  buckling 
stress  of  the  test  specimens  was  estimated  from  an  analysis  for  an  orthotropic  shell 
developed  by  Chao  in  reference  11.  The  analysis  is  similar  to  the  buckling  solution  pre¬ 
sented  by  Cheng  and  Ho  in  reference  12.  Minor  differences  in  the  two  solutions  are  asso¬ 
ciated  with  the  use  of  Timoshenko’s  shell  theory  by  Chao  and  Fltigge's  shell  theory  by 
Cheng  and  Ho.  The  theory  rigorously  accounts  for  the  generally  orthotropic  character  of 
the  tube  helical- layer  configurations.  Structural  stiffnesses  employed  in  the  determina¬ 
tion  of  the  buckling  strength  were  derived  from  the  Halpin-Tsai  equations  and  the  proce¬ 
dures  described  in  appendix  B.  The  buckling  computations  were  performed  with  a  com¬ 
puter  program  based  on  reference  11. 

The  significance  of  boundary  conditions  on  the  predicted  torsional  buckling  stress  of 
the  test  specimens  was  investigated  by  making  buckling  calculations  for  both  clamped  and 
simply  supported  cylinders.  The  boundary  conditions  prescribed  during  buckling  were 
those  usually  employed  in  classical  theory  (for  example,  ref.  13)  and  can  be  written  as 
follows : 


For  simply  supported  cylinders, 
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For  clamped  cylinders, 
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The  results  of  the  computations  are  presented  in  figure  9  where  the  predicted  values 
of  shear  stress  corresponding  to  elastic  torsional  instability  are  compared  with  the  exper¬ 
imental  torsional  strength  of  the  tube.  The  torsional  strength  was  defined  as  either  the 
0.2-percent  offset  yield  stress  Ty  or  the  maximum  shear  stress  rmax  of  speci¬ 
mens  for  which  the  0.2-percent  offset  yield  was  not  defined.  The  buckling  computations 
are  based  on  a  nominal  volume  fraction  and  are  for  a  shell  with  a  length  corresponding  to 
the  unsupported  test-section  length  of  7.6  cm  (3  inches).  The  buckling  results  presented 
correspond  to  0.67  times  the  theoretical  value  obtained  from  Chao's  work.  This  reduction 
in  theoretical  results  has  been  suggested  for  moderately  long  cylinders  in  reference  14, 
and  it  is  believed  to  be  appropriate  for  the  test  specimens.  Although  the  specimens  are 
relatively  thick  walled  shells  (R/t  =  30),  they  appear  to  behave  as  moderately  long  ortho¬ 
tropic  shells,  if  the  curvature  parameter  Z  is  used  as  a  length  measure.  For  ortho¬ 
tropic  shells,  the  parameter  is  a  function  of  structural  stiffness  as  well  as  of  geometry 
(ref.  15)  and  can  be  written  as 


Z  = 


cn;4 

12D22R2 


(4) 


where  is  the  axial  extensional  stiffness  and  D22  is  the  circumferential  bending 

stiffness.  For  the  test  specimens,  Z  ranges  from  70  to  180,  which  is  in  the  moderate- 
length  range.  Studies  of  isotropic  cylindrical  shells  of  similar  Z  (see  refs.  16  and  17) 
have  suggested  that  theoretical  torsion  buckling  loads  should  be  reduced  to  account  for 
imperfection  sensitivity. 

The  results  shown  in  figure  9  indicate  that  strengths  for  tubes  with  only  30°  windings 
and  those  with  only  45°  windings  correlate  with  the  buckling  predictions.  The  observed 
strengths  for  other  helical  winding  angles  and  for  the  tubes  with  both  helical  and  circum¬ 
ferential  windings  are  significantly  lower  than  the  buckling  predictions.  On  the  basis  of 
these  trends,  it  appears  that  most  of  the  specimens  were  not  subject  to  elastic  instability. 
The  strain- reversal  trends  as  indicated  in  table  IV  substantiate  the  buckling  behavior  of 
all  the  tubes  with  only  45°  windings  and  for  one  of  the  tubes  with  only  30°  windings. 

Hence,  it  is  believed  that  tubes  with  only  45°  windings  failed  by  elastic  instability.  It  is 
not  obvious  whether  the  tubes  with  only  30°  windings  buckled  elastically.  The  strain- 
reversal  trends  noted  for  other  tubes  in  tables  IV  and  V  occurred  at  strains  well  beyond 


13 


those  corresponding  to  0.2-percent  offset  yield  stresses  and  are  believed  to  be  the  conse¬ 
quences  of  stiffness  degradation  and  inelastic  buckling  at  loads  corresponding  to  large  non¬ 
linear  strains. 


Yielding  of  Layers 

Stress  analysis.-  The  elastic  stresses  necessary  to  cause  yielding  in  the  layers  of 
the  test  specimens  were  investigated.  By  using  equations  (1)  and  the  orthotropic  shell 
analysis  described  in  appendix  B,  stresses  in  each  layer  were  determined  in  the  cylindri¬ 
cal  surface  coordinates  (x,y  in  fig.  1).  The  stresses  were  then  resolved  into  components 
along  axes  parallel  and  transverse  to  the  fibers  (L,T)  by  employing  the  usual  stress  trans¬ 
formation  equations. 

Results  for  the  average  stresses  ctl>  °T  an(^  tlt  in  each  layer  are  presented 
in  figures  10  and  11.  The  curves  shown  indicate  the  changes  in  these  stresses  with  helical 
winding  angle.  Because  the  membrane  stresses  in  each  layer  are  linear  with  respect  to 
applied  stress  (see  appendix  B),  the  stresses  shown  have  been  normalized  with  respectto 
the  applied  shear  stress  Txy.  In  general,  for  the  torque  shown  in  figure  1  (that  is,  +T 
and  +Txy) >  the  relative  state  of  stress  in  the  helically  wound  layers  at  +oi  and  -a  can 
be  distinguished  in  sketch  (a): 


a. 


>-  ar 


+a  layers 


■a  layers 


Sketch  (a) 


The  circumferential  layers,  of  course,  are  always  in  a  state  of  pure  shear. 

Two  conclusions  can  be  drawn  from  the  results  presented  in  figures  10  and  11.  The 
first  is  that  the  presence  of  transverse  tensile  stress  in  the  -a  layers  is  likely  to  cause 
yielding  in  the  matrix  in  these  layers  before  yielding  of  the  +ot  layers.  The  conclusion 
is  based  on  the  observation  that  tests  of  circumferentially  wound  tubes  in  tension  and 
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compression  have  indicated  that  yielding  occurs  in  tensilely  stressed  tubes  at  significantly 
lower  stresses.  If  yielding  does  occur  in  the  -a  layer,  the  presence  of  compressive 
stresses  in  the  filament  direction  would  be  expected  to  cause  microbuckling  or  wrinkling 
of  the  filaments  since  the  yielded  matrix  presumably  would  have  diminished  capability  to 
support  the  filaments.  On  the  other  hand,  the  +a  layer  is  inherently  more  stable 
because  the  filaments  in  this  layer  are  loaded  in  tension.  Thus,  yielding  of  the  matrix 
because  of  transverse  compression  stresses  is  not  likely  to  precipitate  filament  failure. 
Hence,  it  is  concluded  that  the  -a  layers  are  the  weakest  elements  in  the  composite 
structure. 

The  second  conclusion  is  that  the  critical  stress  interactions  in  individual  layers 
are  associated  with  shear  stresses  and  stresses  transverse  to  the  filament.  From  fig¬ 
ure  10,  it  can  be  seen  that  the  maximum  stress  in  the  filament  direction  ctl  is  at  most 
l|  times  the  applied  shear  stress.  From  table  V  the  maximum  applied  shear  stress 
was  about  200  MN/m2  (30  ksi).  Since  tensile  or  compressive  strengths  of  tubes  with 
axially  oriented  filaments  are  roughly  six  or  seven  times  this  value,  the  influence  of  the 
magnitude  of  stresses  in  the  filament  direction  should  be  relatively  small. 

Orthotropic  yield  criteria.-  With  the  stresses  in  each  layer  of  the  tube  expressed  in 
terms  of  the  stress  ratios  crL//rxy,  uT/?xy,  and  TLT/?xy,  stresses  at  yield  were  deter¬ 
mined  from  various  forms  of  orthotropic  yield  criteria.  The  more  widely  used  ortho¬ 
tropic  yield  criteria  proposed  in  the  literature  can  be  written  in  the  general  form: 


where  Ty  is  the  value  of  the  applied  shear  stress  rxy  at  yielding  of  the  tube,  X  and 
Y  are  yield  strengths  of  a  unidirectional  layer  in  directions  parallel  and  transverse  to  the 
filaments,  and  S  is  the  yield  strength  in  shear  associated  with  the  filament  coordinates 
When  various  values  are  assigned  to  the  coefficients  a  and  b,  equation  (5)  will  reduce 
to  forms  proposed  in  the  literature  as  indicated  in  the  table: 


a 

b 

Source 

Reference 

0 

0 

ANC-18 

18 

1 

0 

Hill 

19 

0 

1 

Norris 

20 

0 

Empirical 

Chamis 

21 

15 


In  the  analysis,  values  assigned  to  the  strengths  X,  Y,  and  S  were  as  follows: 


Strength 

Tension 

Compression 

Shear 

MN/m2 

ksi 

MN/m2 

ksi 

MN/m2 

ksi 

X 

1379 

200 

1379 

200 

Y 

41 

6 

138 

20 

S 

41 

6 

The  strengths  for  Y  are  based  on  several  compression  tests  and  a  single  tension  test  on 
circumferentially  wound  tubes  similar  to  the  test  specimens.  The  strength  S  is  based 
on  the  present  tests  of  tubes  with  only  90°  windings  whereas  X  is  a  conservative  nomi¬ 
nal  value. 

Calculations  for  yield  strengths  indicated  that  terms  involving  X  in  the  yield  cri¬ 
terion  are  negligible  except  for  helical  wrap  angles  near  45°.  This  trend  is  a  consequence 
of  the  stress  distributions  previously  discussed.  Because  of  this  effect,  a  simplified  form 
of  the  yield  criterion 


appeared  to  be  appropriate  for  the  present  study. 

A  comparison  of  the  experimental  torsional  strengths  corresponding  to  the  0.2- 
percent  offset  yield  stress  (or  maximum  stress  when  the  0.2-percent  offset  yield  could 
not  be  defined  (tables  IV  and  V))  and  analytical  yield  stress  predictions  for  a  nominal  vol¬ 
ume  fraction  of  0.63  is  given  in  figure  12.  Data  for  45°  helically  wrapped  tubes  have  been 
omitted  since  buckling  calculations  and  experimental  observations  suggest  that  these  speci¬ 
mens  failed  by  buckling.  The  solid  curves  shown  correspond  to  the  simplified  yield  cri¬ 
terion  (eq.  (6))  for  the  +a,  -a,  and  90°  layers  of  the  cylinders.  The  dashed  curves  cor¬ 

respond  to  the  more  complex  criterion  (eq.  (5))  as  proposed  by  Norris  (a  =  0,  b  =  1)  and 
have  been  included  to  illustrate  the  significance  of  the  X-terms  in  the  yield  criteria.  The 
curves  shown  verify  that  yielding  should  occur  in  the  -a  layer  before  it  occurs  in  the 
+a  layer  because  of  the  weakness  of  the  -a  layer  in  transverse  tension.  Substantial 
differences  exist  between  yielding  predictions  in  the  +a  layers  and  the  -a  and  90° 
layers.  The  experimental  data  would  be  expected  to  lie  between  the  yield  predictions  for 
+a  and  other  layers,  since  local  yielding  in  one  layer  is  likely  to  shift  load  to  other  layers 
and  cause  a  reduction  in  the  yielding  strength  of  the  multilayered  tube. 
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The  data  shown  in  figure  12,  however,  appear  to  correlate  with  the  +o;  curves. 

This  unexpected  result  can  only  be  explained  by  postulating  that  yielding  in  the  -a  layer 
is  a  relatively  mild  phenomenon.  As  discussed  previously,  the  compressive  stresses  in 
the  -a  layer  are  likely  to  cause  microbuckling  or  wrinkling  of  the  filaments.  Except 
for  the  surface  layer,  the  -a  layers  are  sandwiched  between  two  stable  layers  in  which 
at  least  one  layer  (the  +a  layer)  has  filaments  stretched  in  tension.  Presumably,  this 
condition  permits  the  multilayer  composite  to  carry  substantially  larger  loads  after  initial 
local  yielding.  The  large  nonlinear  strains  observed  in  the  helically  and  circumferentially 
wound  specimens  also  demonstrate  that  unusual  behavior  can  occur  in  complex  multi¬ 
layered  configurations.  In  the  present  calculations,  if  yielding  in  the  -a  layers  and  90° 
layers  of  the  tubes  is  ignored,  the  present  data  could  be  correlated  reasonably  well  with 
orthotropic  yield  criteria.  Alternatively,  the  data  could  also  be  correlated  with  analysis 
if  larger  values  of  the  transverse  tensile  yield  stress  Y  and  the  shear  stress  S  were 
employed  in  the  analysis  so  that  yielding  in  the  helical  (+Q!  and  -o')  and  circumferential 
layers  was  coincident.  The  higher  values,  however,  are  believed  to  be  unrealistic  and  are 
not  advisable  for  design.  Thus,  as  applied  herein,  the  orthotropic  yield  criteria  appear  to 
give  unduly  conservative  predictions  of  the  behavior  of  the  test  specimens. 

The  conservativeness  of  similar  applications  of  orthotropic  yield  criteria  has  also 
been  noted  in  correlating  compression  data  (ref.  22)  for  glass-epoxy  tubes.  The  present 
study  and  referenced  studies  raise  serious  questions  as  to  the  engineering  applicability  of 
the  orthotropic  yield  criteria  based  on  elastic  stress  analysis  to  estimate  compressive  and 
shear  strengths  of  multilayered  composites.  It  appears  that  a  more  intimate  knowledge  of 
the  postyielding  behavior  of  composite  structures  is  required  before  reasonable  estimates 
of  maximum  strengths  can  be  made.  Some  recent  studies  (refs.  23  and  24)  have  suggested 
that  orthotropic  yield  criteria  can  be  applied  to  study  progressive  layer  failures.  Most  of 
the  successes  with  postyielding  analyses,  however,  have  been  obtained  with  layers  loaded 
in  tension;  thus,  applications  to  shear  and  compressive  loadings  deserve  further  study. 

CONCLUDING  REMARKS 

The  results  of  torsion  tests  to  failure  conducted  on  36  multilayered  filament-wound 
tubes  have  been  reported.  Shear  stress-strain  data  obtained  during  the  tests  indicated 
various  degrees  of  nonlinearity  (yielding)  depending  on  the  helical  orientation  of  filaments 
in  the  tubes.  Shear  moduli  based  on  the  initial  slope  of  the  stress-strain  curves  were  in 
good  agreement  with  analytical  predictions  from  a  semiempirical  theory  employing  constit 
uent  material  properties.  Fracture  of  helically  wound  tubes  appeared  to  be  accompanied 
by  separation  lines  parallel  to  filaments  stressed  in  tension.  Fracture  patterns  in  tubes 
with  interspersed  helically  and  circumferentially  wound  layers  appeared  to  be  a  band 
(rather  than  a  separation  line)  in  which  delamination  was  evident.  Specimens  with 
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interspersed  circumferential  windings  exhibited  the  capability  of  carrying  much  larger 
postyi elding  strains  than  comparable  specimens  with  only  helical  windings.  Specimens 
with  only  45°  windings  exhibited  strain  reversal  near  failure  and  are  believed  to  have 
buckled.  The  maximum  strength  of  these  tubes  were  in  agreement  with  torsional  buckling 
predictions  based  on  orthotropic  theory  when  theoretical  values  were  reduced  to  account 
for  imperfection  sensitivity.  Shear  stress  data  at  yield  (as  defined  by  0.2-percent  offset 
strain)  or  maximum  load  were  compared  with  yield  strengths  obtained  from  orthotropic 
yield  criteria  and  elastic  stress  analysis.  Computed  results  indicate  that  the  yield  cri¬ 
teria  do  not  adequately  explain  the  behavior  of  helical  layers  in  the  test  specimens  which 
are  loaded  in  tension  transverse  to  the  filaments.  Data  could  be  correlated  only  when 
yielding  of  certain  layers  was  ignored  or  when  unrealistically  large  values  of  transverse 
tensile  and  shear  strengths  were  employed  in  the  strength  analysis. 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Hampton,  Va. ,  May  26,  1971. 


APPENDIX  A 


CONVERSION  OF  U.S.  CUSTOMARY  UNITS  TO  SI  UNITS 


The  International  System  of  Units  (SI)  was  adopted  by  the  Eleventh  General  Confer¬ 
ence  on  Weights  and  Measures  in  1960.  (See  ref.  9).  Conversion  factors  for  the  units 
used  herein  are  given  in  the  following  table: 


Physical  quantity 

U.S.  Customary 
Unit 

Conversion 

factor 

(*) 

SI  Unit 

Length . 

Temperature . 

Force . 

Density . 

Stress,  pressure .  .  . 
Torque  . 

in. 

(op  +  460) 

lbf 

lbm/ft3 
psi  =  lbf /in3 

in -lbf 

0.0254 

5/9 

4.448 

16.02 

6895 

0.1130 

meters  (m) 
kelvin  (K) 
newtons  (N) 

kilograms  per  cubic  meter  (kg/m3) 
newtons  per  square  meter  (N/m2) 
newton -meters  (N-m) 

^Multiply  value  given  in  U.S.  Customary  Unit  by  conversion  factor  to  obtain  equiva¬ 
lent  value  in  SI  Unit. 


Prefixes  to  indicate  multiple  of  units  are  as  follows: 


Prefix 

Multiple 

kilo  (k) 

103 

mega  (M) 

106 

giga  (G) 

109 

centi  (c) 

10"2 

micro  (p) 

10~6 
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APPENDIX  B 


STRESS  DISTRIBUTION  IN  A  GENERALLY  ORTHOTROPIC 
CYLINDRICAL  SHELL 


Stress  distributions  in  the  test  specimens  were  investigated  by  developing  a  solu¬ 
tion  for  axisymmetric  deformation  in  a  multilayered,  orthotropic,  cylindrical  shell  loaded 
in  torsion  and  axial  compression.  The  orthotropicity  in  the  cylinder  layers  was  such  that 
neither  of  the  principal  axes  of  each  layer  was  assumed  to  be  alined  with  the  cylinder  axis. 
The  details  of  the  general  solution  are  rather  lengthy  and  are  not  presented  herein.  A 
general  conclusion  of  the  investigation  was  that  a  state  of  membrane  stress  is  developed 
under  torsional  and  compressive  loadings  over  a  large  area  of  the  tubes  and  that  effects  of 
restraint  at  the  boundaries  of  the  shell  are  very  localized.  The  results  presented  herein 
are  intended  to  document  those  equations  pertinent  to  analytical  estimates  of  the  composite 
wall  stiffnesses  of  the  tube  and  the  stress  distribution  in  each  of  the  layers  of  the  tube. 

In  the  analysis,  the  elastic  constants  for  each  of  the  layers  of  the  tube  were  taken  to 
be  generally  orthotropic  so  that  Hooke’s  law  for  the  ith  layer  can  be  written  as 


where  the  superscript  T  denotes  the  total  strain  in  the  layer  (bending  +  extension).  The 
elastic  constants  in  the  x,y  coordinate  system  (fig.  1)  are  related  to  properties  perpen¬ 
dicular  and  parallel  to  the  fiber  direction  by  the  following  transformation  equations 
(ref.  25) 
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APPENDIX  B  -  Continued 


As  discussed  in  the  text,  the  elastic  constants  associated  with  the  unidirectional  layer  (L 
and  T  coordinates)  were  estimated  from  the  Halpin-Tsai  equations.  The  stresses  in  the 
ith  layer  of  the  shell  can  be  found  from  the  inverse  of  equation  (Bl)  as 


O’, 


xy 


i  «- 


Bn 

Bl2 

b16 

B12 

B22 

B26 

b16 

b26 

b66 

— *i 


:x 


xy 


(B3) 


-■l 


The  shell  was  considered  to  experience  only  axisymmetric  deformation.  By  using 
the  Donnell- von  Karman  assumptions,  the  strains  appearing  in  equations  (B3)  can  be  writ¬ 
ten  in  terms  of  the  tangential  and  normal  displacements  u,  v,  and  w  as 


exT  =  ex  +  zkk  =  u,x  + 1  wjX2  -  zw;XX 


T  _  f  _  w 

y  _ey_R 


(B4) 


yxyT=yxy  =  vx 


To  employ  shell  theory,  it  is  convenient  to  integrate  equations  (B3)  across  the  shell  wall 
and  introduce  the  following  stress  resultants: 


~NX" 

Oil  C12  Ci6  Kn 

ex 

Ny 

c12  c22  c26  K12 

ey 

Nxy 

^16  ^26  ^66  Kjg 

yxy 

Mx 

K11  k12  K16  D11 

.Kxj 

(B5) 


where  Cjk,  Dj^,  and  K^,  the  structural  stiffnesses,  associated  with  extension,  bending, 
and  bending- extension  coupling,  respectively,  are  defined  by 

^Cjk,Kjk,Dj^  =  ^  ^  ^jk^»z,z2)dz  (B6) 

By  use  of  the  Donnell  theory,  the  equilibrium  equations  for  the  shell  can  be  written  as 


N*,x  =  0 


(B7) 
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APPENDIX  B  -  Continued 


Nxy,x  =  0  (B8) 

“MX  XX  +  ~R  "  Nxw,xx  =  p  <B9> 

where  p  denotes  internal  pressure.  The  general  case  in  which  an  axial  compressive 
load  Nx  and  shear  load  are  applied  at  the  boundaries  was  investigated.  For  this 

case  equations  (B7)  and  (B8)  can  be  integrated  so  that 


Nx  =  -Nx 
Nxy  =  Nxy 


(BIO) 

(Bll) 


Equation  (BIO)  together  with  the  first  row  of  matrix  equation  (B5)  imply  that 


K 


ex  ~ 


11 


Cll 


w 


XX 


C12  w  £16  _£x_ 

Cll  R  Cn  >x  Cn 


(B12) 


If  equations  (Bll)  and  (B9)  are  expressed  in  terms  of  displacement,  and  equation  (B12)  is 
used  to  eliminate  ex,  the  following  coupled  equations  in  v  and  w  result: 


with 


Alv,x  +  a2w,xx  +  a3w  ^  P1 

a2v,xxx  +  a3v,x  +  a4w,xxxx  +  A5w,xx  +  A6W  =  p2 


A1  =  c66  - 


C162 


11 


A2=SMU.k16 


a4  =  D11 


kid 

cn 


A5=if^r-KirNx 


11 

\ 

—  —  Ci  c 

P1  =  Nxy  +  Nx  0^ 

nxc12 

Po  =  ■*■■  +  p 

2  RC11 


(B13) 

(B14) 
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APPENDIX  B  -  Continued 


By  use  of  the  usual  techniques  for  solving  simultaneous  ordinary  differential  equations,  a 
general  solution  of  equations  (B13)  and  (B14)  can  be  found  in  terms  of  hyperbolic  and  trigo¬ 
nometric  functions.  If  the  boundary  conditions  are  applied  at  x  =  where  l  is  the 
length  of  the  shell,  the  nontrivial  solution  becomes  antisymmetric  in  v  and  symmetric 
in  w  and  can  be  expressed  as 


v  =  Bi  sinh  /3x  cos  yx  +  B2  cosh  fix  sin  yx  + 


w  =  Cj  sinh  /3x  sin  yx  +  C2  cosh  |3x  cos  yx  + 


P]A6  -  P2A3 

9  X 

(B15) 

AiA6  -  A32 

P2A1  "  P1A3 

O 

(B16) 

a1a6  "As2 

where 


with 


Sj  =  AjA-4  -  A22 

52  =  A1A5  -  2A2A3 

53  =  A]_A6  -  A32 


Calculations  employing  the  general  solutions  (eqs.  (B15)  and  (B16))  with  clamped  or 
simply  supported  boundary  conditions  indicate  that  in  a  multilayered  tube  with  alternating 
layers  at  +a  and  -a  the  boundary  conditions  are  very  localized.  Hence  the  membrane 
displacements  contained  in  the  particular  solutions  of  equations  (B15)  and  (B16)  should  be 
good  approximations  with  which  to  determine  the  state  of  stress  in  regions  of  the  shell 
away  from  the  ends.  For  the  configurations  of  the  present  investigation, 


c16  =  c26  =  0 


Thus  for  the  tube  loaded  only  in  torsion,  the  membrane  displacements  are 
^xy 

v  =  •pr~*~ x 
c66 

w  =  0 


(B17) 


(B18) 
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APPENDIX  B  —  Concluded 


By  using  the  definition  of  the  stress  resultants  (B5),  the  effective  shear  modulus  of 
the  multilayered  wall  is  defined  as 

—  _  ^xy/t  _  ^66 


G 


xy 


(B19) 


xy 


where  t  is  the  total  thickness  of  the  shell  wall.  From  equations  (B18),  the  stresses  in 
any  layer  of  the  cylinder  can  be  computed  by  employing  equations  (B3)  and  (B4). 

For  configurations  similar  to  those  of  the  present  investigation  which  are  loaded 
only  in  axial  compression,  the  membrane  displacements  are 


v  =  0 


w  = 


12 


NxR 


Cll 


■22 


c122 

C11 


(B20) 


J 


(B21) 


The  effective  axial  modulus  Ex  of  the  multilayered  shell  is  given  by 

The  effective  Poisson's  ratio  associated  with  axial  compressive  loading  can  be  written  as 

ey  _  c12 


My  —  ” 
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(B22) 
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TABLE  I.-  TUBE  CONSTITUENT  PROPERTIES 


(a)  International  System  of  Units 


Material 

E 

GN/m2 

kg%i3 

Type  E-HTS  12  and  20  end 
glass  roving 

72.40 

0.23 

2.57 

ERL  2256/MPDA  epoxy  resin 
Cure:  2  hours  at  394  K 

4.37 

0.40 

1.23 

(b)  U.S.  Customary  Units 


Material 

E, 

ksi 

p, 

lbm/cu  in. 

Type  E-HTS  12  and  20  end 

10  500 

0.23 

0.0930 

glass  roving 

ERL  2256/MPDA  epoxy  resin 

635 

0.40 

0.0445 

Cure:  2  hours  at  250°  F 

TABLE  II.-  GEOMETRY  AND  VOLUME  FRACTIONS  OF  HELICALLY  WOUND  TUBES 
[Tube:  inside  diameter,  7.62  cm  (3.00  in.);  test  section  length,  7.62  cm  (3.00  in.)] 


Specimen 

Winding  angle,  a 

deg 

t 

vf 

vm 

vv 

Nominal 

+ 

- 

cm 

in. 

1 

15 

15 

19 

0.160 

0.0631 

0.481 

0.495 

0.024 

2 

15 

14 

18 

.160 

.0631 

.491 

.488 

.021 

3 

15 

13 

17 

.151 

.0593 

.537 

.445 

.018 

4 

15 

20 

17 

.148 

.0584 

.524 

.463 

.013 

5 

30 

32 

30 

.144 

.0567 

.580 

.388 

.032 

6 

30 

30 

31 

.144 

.0569 

.577 

.375 

.048 

7 

30 

31 

32 

.168 

.0660 

.565 

.402 

.033 

8 

30 

32 

30 

.168 

.0663 

.558 

.398 

.044 

9 

45 

48 

46 

.136 

.0534 

.607 

.325 

.068 

10 

45 

46 

48 

.131 

.0516 

.609 

.328 

.063 

11 

45 

48 

45 

.130 

.0511 

.621 

.339 

.040 

12 

45 

47 

50 

.136 

.0534 

.627 

.336 

.037 

13 

60 

61 

62 

*.163 

*.0641 

.651 

.278 

.071 

14 

60 

61 

61 

*.157 

*.0620 

.667 

.275 

.058 

15 

60 

63 

61 

.139 

.0547 

.651 

.296 

.053 

16 

60 

61 

61 

.141 

.0555 

.651 

.295 

.054 

17 

75 

76 

75 

*.101 

*.0396 

.620 

.334 

.046 

18 

75 

77 

77 

.102 

.0402 

.637 

.342 

.021 

19 

75 

77 

77 

*.161 

*.0633 

.585 

.367 

.048 

20 

75 

76 

75 

*.164 

*.0646 

.584 

.374 

.042 

21 

90 

90 

.142 

.0559 

.633 

.324 

.043 

22 

90 

89 

.143 

.0562 

.630 

.312 

.058 

23 

90 

90 

*.128 

*  .0506 

.636 

.297 

.067 

24 

90 

89 

*.124 

*  .0490 

.656 

.273 

.071 

*  Scatter  in  t  >  ±5  percent. 
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(a)  International  System  of  Units 


30 


Significant  strain  reversal  at  failure. 


TABLE  IV.-  TEST  RESULTS  FOR  HELICALLY  WOUND  TUBES  -  Concluded 
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Significant  strain  reversal  at  failure. 


TABLE  V.-  TEST  RESULTS  FOR  HELICALLY  AND  CIRCUMFERENTIALLY  WOUND  TUBES 
(a)  International  System  of  Units 


Nominal  winding  angle, 

Shear  modulus, 

Gxy,  MN/m2 

Experimental  offset 

a,  deg 

Experimental 

Analytical 

yield  stress,  r y,  MN/m2 

25 

15,  90 

9  600 

9  350 

60.0 

26 

15,  90 

9100 

9  887 

54.5 

27 

15,  90 

9  000 

9  584 

60.7 

28 

15,  90 

9  000 

9  522 

62.7 

29 

30,  90 

11  900 

11377 

104.1 

30 

!  30,  90 

11300 

11  515 

104.1 

31 

30,  90 

11200 

11439 

111.0 

32 

30,  90 

11300 

11239 

115.8 

33 

45,  90 

13  400 

12  880 

159.3 

34 

45,  90 

13  200 

12  728 

157.2 

35 

45,  90 

13  800 

13  445 

157.9 

36 

45,  90 

14  200 

13  294 

163.4 

(b)  U.S.  Customary  Units 


Experimental_offset 
Experimental  I  Analytical  I  V^ld  stress,  Ty,  ksi 


"max’ 
in. -kips 


Significant  strain  reversal  at  failure. 


Figure  3.-  Torsion  test  setup. 


i 

- 


/ 

v 


••  • 

■■ . 


TENSION 


Nominal  winding  angles 
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Figure  6.-  Typical  appearance  of  tubes  at  conclusion  of  torsion  tests. 
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(a)  Tubes  with  only  helical  windings,  t  =  0.13  cm  (0.05  in.);  y  =  30. 


a,  deg 


R 

(b)  Tubes  with  helical  and  circumferential  windings,  t  =  0.15  cm  (0.06  in.);  —  =  2! 

Figure  9.-  Comparison  of  experimental  torsional  strengths  of  tubes  with  torsional 
buckling  predictions.  v*  =  0.63;  =  3.0;  ?s  =  3.5;  4  =  2. 


Figure  10.-  Variation  of  stress  with  helical  winding  angle  in  tubes  with  only 
helical  windings,  =  0.63;  ^  =  3.0;  =  3.5. 


(a)  Normal  stresses  in  helically  wound  layers. 

Figure  11.-  Variation  of  stress  with  helical  winding  angle  in  tubes  with  helical 
and  circumferential  windings.  Vj  =  0.63;  £t  =  3.0;  =  3.5. 
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Yield  criterion 

-  Simplified  (eq ,  (6)) 

- Norris  (ref.  20) 


a,  deg 


(a)  Tubes  with  only  helical  windings. 


a,  deg 


(b)  Tubes  with  helical  and  circumferential  windings. 

Figure  12.-  Comparison  of  experimental  torsional  strengths  of  tubes 
with  analytical  predictions  based  on  orthotropic  yield  criteria. 
vf  =  0.63;  =  3.0;  £s  =  3.5. 
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